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We study the next-to-minimal supersymmetric standard model. We consider soft 
supersymmetry breaking parameters, which are induced by the mirage mediation 
mechanism of supersymmetry breaking. We concentrate on the mirage mediation, 
where the so-called mirage scale is the TeV scale. In this scenario, we can realize 
the up- type Higgs soft mass of 0(200) GeV, while other masses such as gaugino 
masses and stop masses are heavy such as 1 TeV or more. Cancellation between the 
effective /U-term and the down-type Higgs soft mass ameliorates the fine-tuning in 
the electroweak symmetry breaking even for u = 0(500) GeV. The mixing between 
the doublet and singlet Higgs bosons is suppressed by (X/k) tan -1 (3. Then the 
lightest doublet Higgs mass naturally reaches 125 GeV lifted by the new quartic 
coupling. The higgsino and singlino are light and their linear combination is the 
lightest superparticle. 



Abstract 



1 Introduction 



Supesymmetric extension is a good candidate for physics beyond the standard model 
(SM). The minimal supersymmetric standard model (MSSM) is the simplest extension. 
The MSSM is quite interesting because of its minimality and its detailed studies have 
been done for several aspects. 

However, the MSSM has the fine-tuning problem. Within the framework of the MSSM, 
the Z-boson mass, mz, is obtained as 



2 

~ ~m 2 Hu - |/x| 2 , (1.1) 
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where m 2 H is the soft supersymmetry (SUSY) breaking scalar mass squared of the up- 
sector Higgs field and /x is the supersymmetric mass. The radiative corrections on m 2 H 
are obtained as m 2 Hu ~ — m| ~ ~M 2 , where mi and M3 denote the stop and gluino 
masses, respectively. In most of cases, the stop and the gluino masses are much larger 
than mz- Thus, we need fine-tuning between m 2 H and |/z| 2 to realize the correct value 
of mz- Furthermore, it is required that mi = 0(1) TeV or larger in order to obtain 
the Higgs mass such as m^ w 125 GeV which is recently reported by ATLAS and CMS 
collaborations (21 E]- 

The mirage mediation is one of the interesting mediation mechanisms of SUSY break- 
ing [H |5j [6] . The mirage mediation is a mixture of the modulus mediation |7j and the 
anomaly mediation [8] with a certain ratio. In particular, it was pointed out that the 
TeV-scale mirage mediation can ameliorate the above fine-tuning problem of the MSSM 
[HI Ell [EE]- In the TeV-scale mirage mediation, the above radiative corrections on m 2 H and 
the anomaly mediation contributions are canceled each other. Then, the value of \m% | at 
the electroweak scale can be smaller than stop and gluino masses. The TeV scale mirage 
mediation also leads several phenomenologically interesting aspects [12] because its SUSY 
particle spectrum is quite compressed. 

The next-to-minimal supersymmetric standard model (NMSSM) is the extension of 
the MSSM by adding a singlet S [13] (see for review e.g. [E]). Here, we also impose 
the Z 3 symmetry. The NMSSM does not have the //-term, {iH u Hd, in the superpotential, 
where H u and Hd denote the up and the down-sector Higgs superfields, respectively. On 
the other hand, the term \SH u Hd is allowed in the NMSSM superpotential. After S 
develops its vacuum expectation value (VEV), the effective /i-term is generated. That 
gives us a solution for the so-called //-problem [15J. The NMSSM is also interesting in the 
light of the recent indication of the relatively heavy Higgs boson reported by ATLAS and 
CMS, endowed with an additional Higgs self-coupling. In the Higgs sector, the doublet 
Higgs and singlet fields mix each other. Then, the Higgs sector in the NMSSM has a quite 
rich structure. 



The NMSSM also leads the same relation as (1.1 ), and the NMSSM has the fine-tuning 
problem similar to the one in the MSSM. Thus, it is interesting to apply the TeV-scale 
mirage mediation to the NMSSM. In this paper, we study the NMSSM with the soft 
SUSY breaking terms induced through the TeV-scale mirage mediation pp. 
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This paper is organized as follows. In section 2, we give a brief review on the mirage 
mediation, and the TeV scale mirage mediation. In section 3, we apply the TeV scale mi- 
rage mediation to the NMSSM, and study its spectrum. Section 4 is devoted to conclusion 
and discussion. In Appendix A, we show explicitly initial conditions of soft parameters, 
which are induced through the mirage mediation in the NMSSM. 



2 TeV- scale mirage mediation 

Here, we review briefly the mirage mediation [I]. The mirage mediation is the mixture 
between the modulus mediation and the anomaly mediation. Then, the gaugino masses 
are obtained as 

M a = M + ^&^, (2.1) 

where g a and b a are the gauge couplings and their (3 function coefficients, and m 3 / 2 denotes 
the gravitino mass. The first and the second terms in the right-hand side of Eq. (2.1) 
correspond to the contributions due to the modulus mediation and the anomaly mediation, 
respectively. 

Similarly, we obtain the so-called A-terms corresponding to the Yukawa couplings, 
i/ijk, and the soft scalar masses m ; as 



Aijk(M GUT ) = a ijk M Q - (7» + 7j + 7*.) 



m 



3/2 



87T 
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m*(M GUT ) = Q M 2 -7,(^) 2 -^M ^, (2.2) 



where 



li 



2 

jk 



a jk 



7* = (2-3) 

Here, Cv? denotes the quadratic Casimir corresponding to the representation of the 
matter field <ft l . In addition, and q parametrize the A-term and the scalar mass 
squared generated through the modulus mediation in the unit of the universal gaugino 
mass, M . These coefficients are determined by modulus-dependence of the Kahler metric 
as well as the Yukawa coupling. One can write 

a = cf ree) + 5c r (loop) . (2.4) 

Here, cf me ^ is calculated from the tree-level Kahler metric of the matter field 0* and they 
are ratios of small integers including and 1 jTJ [161 HI]- 111 addition, 5cf oop ' ) is obtained 
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with the one-loop Kahler metric of the matter field, but such a loop correction to the 
Kahler metric depends on the detail of the ultraviolet-model and is hard to calculate (see 
e.g. Ref. [I?]). That is, q is ambiguous at the one- loop level, although such ambiguity is 
subdominant and less important in most of cases. Here, we consider the case with 

aijk = Ci + Cj + c fc . (2.5) 

We input the values of q at Mqut = 2 x 10 16 GeV. 
It is convenient to use the following parameter [5], 

_ ^3/2 , s 

a ~ M ln(M pl /m 3/2 y 1 ' ] 

to represent the ratio of the anomaly mediation to the modulus mediation. Here M p i is 
the reduced Planck scale. 

One of the interesting aspects in the mirage mediation is that the above spectrum 



(2.1) and (2.2) has a special energy scale, that is, the mirage scale, 



M ■ = Mgut (2 7) 

At this scale, the gaugino masses are obtained as [5], 

M a (M mir ) = M . (2.8) 

That is, the anomaly mediation contribution and the radiative corrections cancel each 
other, and the pure modulus mediation appears at the mirage scale. Furthermore, the 
A-terms and the scalar masses squared also satisfy^ 

%fc(M mir ) = (a + Cj + c k )M , m'(M mir ) = qM 2 , (2.9) 

if the corresponding Yukawa couplings are small enough or if the following conditions are 
satisfied, 

a ijk = Ci + cj + c k = I, (2-10) 
for non- vanishing Yukawa couplings, [5] . 



When a = 2, the mirage scale M mir is around 1 TeV. Then, the above spectrum (2.8) 



and (2.9) is obtained at the TeV scale. That is the TeV scale mirage mediation scenario. 
In particular, there would appear a large gap between Mo and the scalar mass rrii with 
Cj ~ 0. We will apply the TeV scale mirage scenario to the NMSSM in the next section. 

In the TeV scale mirage scenario, the stop mass squared becomes negative at high 
energy [T5], while it is positive at low energy below 10 6 GeV. Thus, the vacuum which 

1 The scaler masses at the Mirage scale can be modified due to the U(l) tadpole contribution in the 
renormalization group running when the different values of ch u and cu d are chosen. However, such a 
modification is small and can be included in ambiguities of Ci if couplings are small. See [5] for detailed 
discussions. We include this contribution in our numerical analysis. 
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breaks the electroweak symmetry at the electroweak scale might be a local minimum, but 
instead there would be a color and/or charge breaking vacuum with field values larger 
than 10 6 GeV. Here, we assume the thermal history of the Universe such that field values 
remain around the origin until the temperature reaches the electroweak scale. In addition, 
we need to confirm that the tunnelling rate is small enough, i.e. less than the Hubble 
expansion rate. In Refs. [Til], it has been shown that such a rate is small enough, as long 
as the squark/slepton masses squared are vanishing or positive around 10 4 GeV. This 
condition is satisfied in our TeV scale mirage mediation scenario. 

3 TeV scale mirage in NMSSM 

In this section, we apply the TeV scale mirage mediation scenario to the NMSSM. 
3.1 NMSSM 

Here, we briefly review on the NMSSM, in particular its Higgs sector before we apply 
the TeV scale mirage mediation scenario to the NMSSM. In the NMSSM, we extend the 
MSSM by adding a singlet chiral multiplet S and imposing a Z3 symmetry. Then, the 
superpotential of the Higgs sector is written as 

W mggs = -XSH u H d + ^S 3 . (3.1) 

Here and hereafter, for S, H u and H d we use the convention that the superfield and its 
lowest component are denoted by the same letter. The full superpotential also includes 
the Yukawa coupling terms between the matter fields and the Higgs fields, which are the 
same as those in the MSSM. 

The following soft SUSY breaking terms are induced in the Higgs sector, 

Koft = m 2 Hu \H U \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 - XA x SH u H d + ^A K S 3 + h.c. (3.2) 

Then, the scalar potential of the neutral Higgs fields is given as 

V = X 2 \S\ 2 (\H d \ 2 + \H u \ 2 ) + \KS 2 -XH d H° u \ 2 + V D 

+m 2 H jH u \ 2 + m 2 H \H d \ 2 + m 2 s \S\ 2 - XA x SH u H d + ^A K S 3 + h.c, (3.3) 

with 

vb = ^+^)(iflsr-w 2 ) 2 , (3.4) 

where gi and g 2 denote the gauge couplings of U(l)y and SU(2). 
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The minimum of the potential is obtained by analyzing the stationary conditions of 
the Higgs potential, 

d \^ 1 

\ 2 v cos (3(s 2 + v 2 sin 2 (3) — Xkvs 2 sin (3 H — g 2 v 3 cos (3 cos 2(3 



dH° d ^~ ' I 

,2 



+ rn Hd v cos/3 — \A x vs sin/3 = 0, (3.5a) 
= A 2 t> sin /3(s 2 + w 2 cos 2 /3) — Xkvs 2 cos /3 — -rg 2 v 3 sin /3 cos 2/3 

u 

+ m 2 Hu vsm(3 — \A x vs cos (3 = 0, (3.5b) 



91/ 

where g 2 = g 2 + #2- Here, we denote VEVs as 



\ 2 sv 2 + 2k 2 s 3 - \KV 2 ssm2f3 + m 2 s s AA A t> 2 sin 2/3 + kA k s 2 = 0, (3.5c) 



^ = (I^T> + (I#T>, tan/3 = s = <5>. (3.6) 

Using the above stationary conditions, we obtain the Z boson mass m| = \g 2 v 2 as 

2 l-cos2/3 2 l + cos2/3 2 2 
m z = —7^ m H M —7^7, m H d ~ ^ , (3-7) 



where /i = As. For tan/3 ^> 1, this equation becomes 

2 

m| ~ -2m 2 H + - 2/i 2 . (3.8) 

z Hu tan 2 /3 H * ^ y ' 

This relation is the same as the one in the MSSM. Indeed, when we neglect the second 



term in the right-hand side, the above relation is nothing but Eq. (1.1). Thus, the 
natural values of and would be of (9(100) GeV. Furthermore, the natural value 

of \rriH d \/ tan/3 would be of (9(100) GeV or smaller. Alternatively, and |m#J/tan/3 
could be larger than (9(100) GeV when fi 2 and m 2 H J tan 2 f3 are canceled each other in the 
above relation at a certain level. Even in such a case, I^hJ would be naturally of (9(100) 
GeV. On the other hand, other sfermion masses as well as gaugino masses must be heavy 
as the recent LHC results suggested. To realize such a spectrum, we apply the TeV scale 
mirage mediation in the next section, where we take ch u = to realize a suppressed value 
of \m,H u \ compared with M . 

3.2 TeV scale mirage mediation in NMSSM 

Here, we study the TeV scale mirage mediation scenario in the NMSSM. Soft SUSY 



breaking terms are obtained through the generic formulas (2.1) and (2.2) with taking 
a = 2. For concreteness, we give explicit results of all the soft SUSY breaking terms for 
the NMSSM in Appendix |A| We concentrate on the Higgs sector as well as gauginos and 
stops. 
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We consider the following values of q, 



(tree) _ , (tree) _ n (tree) _ „ (tree) _ (tree) _ 1 /q q\ 

C H d - L , C // u - U ' C 5 - U ' C *l - C t R - 2' ^- y J 

for H d , H u , S, and left and right-handed (s)top fields, respectively. This is the same 

assignment as the pattern II in Ref. [UJ for the MSSM except for c$- Then, the soft 
parameters due to only modulus mediation contribution are given by 

\Af )modulus ( y 4.A)modulus -^0; (-^-k) modulus 0, 

( m # d )modulus = Mq, (W^) mo dulus = ( m f fl ) modulus = n^O' (3.10) 
( m H u ) modulus = (^s) modulus = 0, 



when we neglect 5cf oop \ The above assignment of Cj (|3.9l) satisfies the condition, (12.101) 



for the top Yukawa coupling and the coupling A, but not for the coupling k. However, we 
do not consider a large value of k to avoid the blow-up of k and A as will be shown later. 
Thus, we obtain the following values, 

A t ^A x ^ M , (3.11) 
m\ d « Ml ml « m} R « ^M 2 , 

up to (9(k 2 /87t 2 ) at the TeV scale. Note that 5cf oop ^ has negligible effects for these values. 

Similarly, we can obtain the values of A K , |m_ej and \m s \ at the TeV scale, however 
those are suppressed compared with M . For such suppressed values, sub-leading correc- 
tions e.g. the one-loop correction on the Kahler metric are not negligible anymore. That 
introduces one-loop ambiguity into the model. Including such corrections, at the TeV 
scale we obtain 

m 2 Hu - 5c£ op) M 2 , ml « 6c^Ml (3.12) 

with 5c£ op) ,5c5j oop) = C(1/8tt 2 ). Note that similar to Eq.|3~ll|, Eq.(|3.12| also includes 
corrections of 0(k 2 Mq /87r 2 ) due to the violation of the mirage unification by k. That is, 
we obtain m 2 Hu = 0,m| = up to C(M 2 /8tt 2 ) at the TeV scale. Similarly, at the TeV 
scale we can obtain, 

A K = 0, (3.13) 

up to O(M /8tc 2 ). Because of such ambiguity, we use A K as a free parameter, which 
must be small compared with M . In addition, we determine the values of Tn 2 H , m 2 s and 



/i (= As) at the electroweak scale from the stationary conditions, (3.5), where we use the 
experimental value mz = ~^9 V — 91.19 GeV and tan/3 as a free parameter. 

Through the above procedure, the parameters, m 2 H , m 2 s and /i, at the electroweak 
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scale are expressed by tan/3, m 2 Hd , A x as follows, 

P W 2(1- Stan J?) I J 

ml = -2 



2 tan 2 g - 1 / m\ d 2 _ m| 

H " tan 2 /3 Vtan 2 /3-l ^ ' 2 



where, 



m\ d (1 - f tan /3) J tan 2 /2A 2 tan 2 /3 - 1 



X ^tan 2 /? \ 1 + tan 2 /3 + l V (? 2 2 tan 2 J m\) ' ( ' U '' j) 

For tan/3 ^> max(l, k/X), these parameters are approximated as, 



m | „ _ 2 (*)' ( -!<-) - («) .4„ f ^SL-A + 2^4m|, (3.16b) 
5 VA/ V^A tan /V ^ V^a tan/3 y g 2 m 2 Hd z ' v ; 

mH ^^-^t^~r^- (36) 

When tan/3 = (9(10), the values of /i, and \ms\ are smaller than Mq by the factor 

tan/3 because m# d ~ — Mq. Thus, the values of fi and could be of (9(100) 

GeV while the other masses of the superpartners are of 0(M O ) = 0(1) TeV. Then, the 
fine-tuning problem can be ameliorated. Furthermore, one can see that the first and 
the second terms in the last equation cancel each other for our choice of q. The next 
leading contributions are of 0{m 2 Hd / tan 4 /3) or 0{m 2 Hd ji/ tan 2 (3A\). Thus, m\ is almost 
determined by m 2 H alone and insensitive to the value of /i. This means that actually 
tan/3 PS 3 is enough to obtain the fine-tuning of | d In m 2 z /d In m 2 H | _1 = m 2 z /2m 2 H = 
0(100)% for Mq ~ 1 TeV. In this case, fi can be as heavy as 0(400) GeV without 
deteriorating the fine-tuning. The origin of this cancellation is easily understood by 
examining the doublet mass matrix, 

C M = - (H d , H* u ) M 2 H (^f^, (3.17) 

where 

M%=( m \ +li2 7 AX \)~( M \V 2 ~ M r). (3.18) 

The modulus mediated contribution Mq cancels in the determinant of the mass matrix, 
det(M 2 H ) ~ //. The heavy mode has mass of 0(M O ), then the mass of the light mode 
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is suppressed as fi 2 /M « \ij tan/3 and a flat direction appears along H u /H d w M //i w 
tan/3. This mechanism was previously observed in [TT] in the context of the MSSM. In 
the NMSSM, the relation rriH d ~ ^4a is well controlled up to the leading contribution of 
the modulus mediation, in contrast to the 5-term (in place of A\) in the MSSM, which 
is a remnant of the fine-tuned cancellation between the terms of 0(777.3/2) and subject to 
uncontrolled corrections. 

In the following section, we show numerically the spectrum of our model. 



3.3 Spectrum 

Here, we study numerically the spectrum of our model. Before showing numerical results, 
we recall our parameters. In our analysis, the free parameters are A, k, tan/3 and A K 
given at the SUSY scale and M given at Mqut- As the SUSY scale, we choose Msusy = 
y / m^m^ ~ Mo/a/2. Using them, we determine all of the soft SUSY breaking parameters 
except for rn 2 Hu , m| and /i, which are determined by using the stationary condition of 
the Higgs potential. Note that in the following numerical analysis we include corrections 
from the 1-loop effective potential, 



M 2 1 ? 

1V1 SUSY I 



(3.19) 



in the stationary conditions (3.5), where Ai represents the mass matrix of our model and 
Str denotes the supertrace. In all of the following numerical analysis, we use A K = —100 
GeV as a typical value of A K . When A and/or k are large at the electroweak scale, they 
blow up below the GUT scale. Thus, we have constrains on large values of A and k by 
requiring that those do not blow up below the GUT scale. 

Figure 1 shows the lightest CP-even Higgs mass rrihi, soft scalar masses of H u and S, 
rriH u and mg, and \i for M = 1200 GeV and tan/3 = 3, in panels (a), (c), (d) and (e), 
respectively. The panel (b) in the figure shows the coupling squared between the lightest 
CP-even Higgs and the Z bosons, g Z zhii as the ra tio to the one in the SM, i.e. 9zzhi/ 9sm- 
The second lightest CP-even Higgs mass is also plotted in panel (f). 

In the figure, the red curve corresponds to the values of A and k at the electroweak 
scale, which blow up at the GUT scale. Thus, we exclude the outside of this curve. The 
gray region around k = corresponds to the region where the tachyonic mode appears in 
the Higgs sector. The yellow region is excluded because the Higgs potential has the false 
vacua studied in Ref. [20], deeper than the realistic vacuum. From these constraints, the 
region with small k/X is disfavored. The gray region around the red curve indicates the 
region where the tree level Higgs mass becomes tachyonic and the iterative procedure we 
employed does not work to estimate the stationary conditions. The quantum corrections 



(3.19) could lift the tachyonic mass, however, we do not calculate it because the region 
has already been excluded by the LEPII bound {m hl > 114.4 GeV). 

The value of fi is around 200—400 GeV, which is consistent with the rough estimation in 



Eq. (3.16a), i.e. \i ~ m 2 H /( A\ tan f3>) ~ M / tan/3. Obviously the expansion in Eq. (3.16a) 



becomes worse for k/X > tan/3, while the expansion holds well for k/X < 0.3 where the 
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value of Im/rJ is around 100 — 200 GeV. In this region the fine-tuning of the parameters for 
the electroweak symmetry breaking is of (9(10)%, even though most of the superpartners 
have heavy masses of O(M ). A large value of /i ~ 400 GeV potentially degrades the 
fine-tuning, however the cancellation renders /i irrelevant to rriz at the leading order of 
the expansion by 1/ tan/3. 



The value of \mg\ is roughly estimated as \ms\ ~ (V2k/X)h in Eq. ( |3.16b| ). Thus, it 
is found that \mg\ ~ 400 GeV for k ~ A and \ms\ increases (decreases) as k/X increases 
(decreases). The value of \ms\ is expected to be suppressed in the TeV scale mirage 
scenario. However, a large value of k/X leads large \ms\ through the stationary condition 
like \ms\ ~ 500 GeV. Such a large value would not be realized in our TeV scale mirage 
mediation scenario, because m| must be suppressed compared with Mq. Thus, the region 
with large k/X and \ms\ > O(M /y/8ir 2 ) is disfavored. In the panel (d), the region with 
\ms\ > M0/V871 2 is filled in pink. It is interesting that the region favored by the TeV scale 
mirage mediation exactly corresponds to the region where the fine-tuning is ameliorated. 



Note that the condition (2.10) also holds well for this region. 

In Fig. l.(a), the lightest CP-even Higgs boson, which dominantly consists of the dou- 
blet scalar here, has a mass around 80-125 GeV. We estimated using NMHDE- 
CAY in the NMSSMTools package [2TJ. We calculated the minimum of the effective 
potential renormalized at Msusy by the iteration starting from the tree-level minimum. 
Then we used the resultant \i as an input of the NMSSMTools. The pole mass m t = 172.9 
GeV was used in our calculation. The qualitative behavior of the mass of the SM-like 
Higgs boson is given by 

A 2 

ml M ~ ml cos 2 2/3 + X V sin 2 2/3 -v 2 (A - Ksin2/3) 2 

K 

+ ^H^) + A{ l -^k))' (3 - 20) 

for ks ^> \A K \, \A\\ in [Hj. The first term and the fourth term are the tree-level contribu- 
tion and the radiative correction in the MSSM, respectively. The second term comes from 
the new quartic Higgs couplings in the NMSSM, and the third term comes from the mixing 
of the doublet scalars with the singlet. The third term is always negative because the mix- 
ing reduces the lightest eigenvalue of the mass matrix. Although the above approximation 
does not apply in the figure, the behavior of the higgs mass can be understood in terms of 
the mixing. The effect of the mixing undermines that of the additional quartic coupling 
A 2 sin 2 2/3 and suppresses the Higgs mass except for the narrow region 4 < X/k < 8. The 
LHC observation, m hl rs 125 GeV, is satisfied only around (A, k) = (0.7,0.11) where the 
mixing vanishes. Such a region may be realized as a quasi-infrared fixed point if A has 
a strong dynamics origin at the GUT scale, while k is suppressed due to the approxi- 
mate Peccei-Quinn symmetry. It is important to stress again that this region is favored 
by the TeV scale mirage mediation and the fine-tuning. The mixing of the doublet and 
the singlet Higgs scalars in the Lagrangian is obtained by rotating the doublet Higgs by 
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/3 (= tan.- 1 v u /v d ) 

AC = -2vX \ji - (A\ + 2ks) cos (5 sin 0] hAS r , (3.21) 

where h = Re(H d ) cos (3 + Re(H u ) sin (3 — v and AS r = Re(S') — s denote dynamical degree 
of freedom of the corresponding Higgs fields. Since our TeV scale mirage mediation 
scenario leads n Mq/ tan/3 and A\ « Mo, the above term is approximated as, 



AC = 4v ^ 



tan/3 



K + o' 1 



A Vtan 2 /3 



/iAS r . (3.22) 



Thus the mixing is automatically suppressed by (k/A) / tan /3 in our scenario. 

The coupling of the lightest CP-even Higgs boson to the Z boson, gzzh^-, is almost 
the same as the one in the standard model in most of the parameter space as expected. 
We show the ratio gzz^ldsu * n the figure, where #sm denotes the Higgs coupling to 
the Z boson in the standard model. The mixing between the doublet and the singlet is 
minimized around ~ 125 GeV, where gzzhx also approaches to its SM value. 

Table 1 shows examples of spectra and gzz^/dsM f° r K ) = (0-10, 0.40), (0.40, 0.10) 
and (0.70,0.11). In the table, for % — 1,2,3 and m ai for i = 1,2 denote three CP- 
even Higgs masses and two CP-odd Higgs masses, respectively. Also, 2 denote two 
eigenvalues of stop masses. Other squark and slepton masses depend on their values of 
Cj, and they are of O(M ) unless Cj = 0. The second lightest CP-even Higgs boson h 2 is 
lighter for (A, «) = (0.40,0.10) and (0.70,0.11) than the one for (A, k) = (0.10,0.40). Its 
dominant component is the singlet Higgs boson S. Its mass decreases as k/X decreases. 
Then, the above behavior of occurs. On the other hand, the dominant component 
of the heaviest CP-even Higgs boson is the down-type Higgs boson H d . Its mass is 
heavy and almost equal to Mo, independent of A and k. By the same reason, the mass 
of the heavier CP-odd Higgs boson a 2 is almost the same as as well as Mo. The 
lightest CP-odd Higgs boson a% is lighter for small k due to the approximate Peccei- 
Quinn symmetry. Also, three gaugino masses are almost the same as M . It might be 
challenging but interesting subject to observe these light extra-Higgs bosons through the 
small mixing with the doublets in LHC and ILC. 

Figure 2 shows the same as Fig. 1 except for Mq = 1500 GeV. Most of mass parameters 
become larger than those for M = 1200 GeV. The lightest CP-even Higgs mass becomes 
heavy due to the heavier stop and the portion of the region with ~ 125 GeV increases. 
Table 2 is the same as Table 1 except M = 1500 GeV. Obviously, all of masses become 
heavier than those in Table 1. The behavior of rrih 2 , and m a2 is the same as the one 
in Table 1. For completeness, we also plot the figure for Mo = 1700 GeV in Fig. 3 and 
list the spectra and the coupling for M = 1700 GeV in Table. 3. 

Figure 4 shows the same as Fig. 1 except for tan — 5. The approximation in 



Eq. (3.16b) and the cancellation work well for k/X < 5. The up-type Higgs mass |m#J 
is around 100 — 250 GeV in most of the parameter space, while the region, k/X < 0.6, 
satisfies \ms\ < 100 GeV and is favored by the TeV scale mirage mediation. In this region, 
the singlet becomes lighter than the doublet and g 2 Z zhil ' 9sm can decrease to O(10) % near 
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(A,«) 


(0.10,0.40) 


(0.40,0.10) 


(0.70,0.11) 


m hl 


105 GeV 


107 GeV 


126 GeV 


m h2 


1261 GeV 


182 GeV 


138 GeV 


m h3 


1700 GeV 


1312 GeV 


1320 GeV 


m al 


512 GeV 


181 GeV 


158 GeV 


m a2 


1260 GeV 


1311 GeV 


1321 GeV 


9zzh-J 9sm 


1.00 


0.95 


0.94 


m Hd 


1.39 x 10 6 GeV" 


1.39 x 10 6 GeV 2 


1.40 x 10 6 GeV 2 


m Hu 


1.29 x 10 5 GeV 2 


5.58 x 10 4 GeV 2 


3.06 x 10 4 GeV" 


m\ 


-1.42 x 10 6 GeV 


-1.02 x 10 4 GeV 


-4.93 x 10 3 GeV 2 


m h 


823 GeV 


823 GeV 


823 GeV 


m h 


849 GeV 


849 GeV 


849 GeV 


n 


217 GeV 


396 GeV 


407 GeV 



Table 1: Spectra for (A, k) = (0.10,0.40), (0.4,0.10) and (0.7,0.11) with tan/3 = 3 and 
M = 1200GeV. 



(A,«) 


(0.10,0.40) 


(0.40,0.10) 


(0.70,0.11) 


m h i 


107 GeV 


110 GeV 


128 GeV 


m h2 


1574 GeV 


231 GeV 


172 GeV 


m h3 


2134 GeV 


1637 GeV 


1641 GeV 


m a i 


572 GeV 


200 GeV 


172 GeV 


m a2 


1573 GeV 


1636 GeV 


1645 GeV 


9zzh-J 9sm 


1.00 


0.98 


0.99 


m Hd 


2.25 x 10 6 GeV^ 


2.18 x 10 6 GeV 2 


2.18 x 10 6 GeV 2 


m Hu 


2.03 x 10 5 GeV^ 


5.02 x 10 4 GeV 2 


3.84 x 10 4 GeV 2 


m\ 


-2.25 x 10 6 GeV 2 


-1.92 x 10 4 GeV 2 


-9.38 x 10 3 GeV 2 


m h 


1023 GeV 


1023 GeV 


1023 GeV 


m h 


1055 GeV 


1055 GeV 


1055 


n 


272 GeV 


495 GeV 


508 GeV 



Table 2: Spectra for (A, k) = (0.10,0.40), (0.40,0.10) and (0.70,0.11) with tan/3 = 3 and 
M = 1500GeV. 
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(X,k) 


(0.10,0.40) 


(0.40,0.10) 


(0.70,0.11) 


m hl 


108 GeV 


112 GeV 


128 GeV 


m h2 


1782 GeV 


264 GeV 


196 GeV 


m h 3 


2422 GeV 


1853 GeV 


1861 GeV 


m al 


609 GeV 


202 GeV 


181 GeV 


m a2 


1782 GeV 


1854 GeV 


1861 GeV 




1.00 


0.99 


1.00 


m Hd 


2.79 x 10 6 GeV" 


2.80 x 10 6 GeV^ 


2.80 x 10 6 GeV^ 




2.62 x 10 5 GeV^ 


6.56 x 10 4 GeV" 


5.09 x 10 4 GeV" 


m s 


-2.91 x 10 6 GeV^ 


-2.66 x 10 4 GeV 2 


-1.31 x 10 4 GeV 2 




1157 GeV 


1157 GeV 


1157 GeV 


m t 2 


1193 GeV 


1193 GeV 


1193 GeV 


n 


308 GeV 


560 GeV 


575 GeV 



Table 3: Spectra for (A, k) = (0.10,0.40), (0.40,0.10) and (0.70,0.11) with tan/3 = 3 and 
Mo = 1700GeV. 



the border of the excluded region by the false vacuum. In such a region the lightest CP- 
even Higgs boson (m^ ~ 80 — 90 GeV) could escape the LEPII bound and the second 
lightest CP-even Higgs boson rrih 2 ~ 125 GeV gives the signal observed in LHC [22] . Note 
that the small mixing with the singlet scalar enhances the second lightest Higgs mass in 
contrast to the lightest Higgs case and makes it easier to obtain the LHC value. The 
production cross section of h 2 ( mainly composed of H u ) through the gluon fusion or the 
vector boson fusion is reduced by O(10) % relative to the SM due to the mixing with 
the singlet. The branching ratio BR{h 2 — > ZZ, WW, 77) is also reduced. This is because 
the width T{h 2 — > bb) does not change due to the suppressed mixing between the heavy 
H d and the singlet, while the widths {h 2 — > ZZ, WW, 77) are reduced due to the mixing 
between H u and the singlet. Such a reduction will be confirmed or excluded by the future 
measurement of the Higgs coupling in LHC and ILC. Table 4 shows the examples of the 
mass spectra and g 2 ZZh J glu for (A, k) = (0.10,0.40), (0.40,0.10) and (0.70,0.11). 

When we increase tan/3, \fi\ decreases as expected from the rough estimation, |/i| ~ 
Mq/ tan/3 in Eq. (3.16a). For example, the value tan/3 = 10 leads \i = 100 GeV or less 
for M ~ 1 TeV, which is excluded by the chargino mass bound by LEPII. The value of 
\m s \ also decreases as tan/3 increases. In opposite view, this means that tuning of \m s \ 
(or 5cf oop ^) increases to obtain large tan /3 for fixed (A, k) and small tan j3 is favored in 
our scenario. 

Finally we comment on the fermionic sector. The singlino mass is estimated as 2k[x/\. 
Since we have /1 around 200 — 500 GeV, the singlino can also be light. Note that the region 
with small k/A is excluded by the appearance of the tachyonic modes and/or the false 
vacua. Thus, we can not lead the singlino much lighter than the higgsino. Also, large 
k/X leads large \ms\, which could not be derived in our TeV scale mirage scenario. Thus, 
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(A,«) 


(0.10,0.40) 


(0.40,0.10) 


(0.70,0.11) 


m hl 


113 GeV 


79 GeV 


65 GeV 


m h 2 


1136 GeV 


126 GeV 


130 GeV 


m h3 


1209 GeV 


1222 GeV 


1228 GeV 


m al 


420 GeV 


138 GeV 


121 GeV 


m a2 


1205 GeV 


1221 GeV 


1227 GeV 


9zzhi/9su 


1.00 


0.24 


0.13 


m Hd 


1.39 x 10 6 GeV 2 


1.39 x 10 6 GeV 2 


1.39 x 10 6 GeV 2 


m Hu 


4.91 x 10 4 GeV 2 


2.00 x 10 4 GeV 2 


1.97 x 10 4 GeV 2 


m 2 s 


-6.23 x 10 5 GeV 2 


8.11 x 10 2 GeV 2 


4.85 x 10 3 GeV 


mi, 


822 GeV 


822 GeV 


822 GeV 


m h 


847 GeV 


847 GeV 


847 GeV 


/i 


146 GeV 


228 GeV 


232 GeV 



Table 4: Spectra for (A, k) = (0.10,0.40), (0.40,0.10) and (0.70,0.11) with tan/3 = 5 and 
Mo = 1200GeV. 



singlino much heavier than the higgsino is disfavored. Thus, both masses of the higgsino 
and singlino would be of the same order. Since all of gaugino masses as well as squark and 
slepton masses are much heavier, the lightest superparticle would be a linear combination 
between the higgsino and singlino, depending on k/X. The gravitino is quite heavy such 
as m 3 / 2 ~ 47r 2 M , as already known in the mirage mediation mechanism [H El E]- 

The lightest neutralino in this model behaves like the linear combination between the 
higgsino and bino in the MSSM and could reproduce the observed abundance of the cold 
dark matter assuming the thermal relic saturates it [23]. While the extra Higgs bosons play 
an important role in the direct detection of them, which could be significantly different 
from the results in the MSSM [21]. The detailed study of the phenomenology of the 
TeV scale mirage mediation in the NMSSM is interesting for its distinct mass spectrum, 
however, beyond the scope of this work [25] . 

4 Conclusion 

We have studied the NMSSM with the TeV scale mirage mediation. The region with 
large k/X requires a large value of \ms\ to satisfy the stationary conditions of the Higgs 
potential. Such a large value could not be realized in our TeV scale mirage scenario and 
therefore such parameter region is disfavored. In the favored region, it is found that we 
can realize |m#J ~ 200 GeV, while other masses are heavy such as 1 TeV. Then, the 
fine-tuning problem is ameliorated. The cancellation between the effective /i-term and 
the down-type Higgs soft mass reduces the sensitivity of fi to mz and \i = (9(500) GeV 
is possible without significantly deteriorating the fine-tuning. The mixing between the 
light doublet and the singlet is suppressed by (k/X) tan -1 (5. For small tan /3 the lightest 
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CP-even Higgs is mainly the doublet and its mass reaches 125 GeV without suppression 
by the mixing. When we increase M , the Higgs mass slowly increases, and also the 
value of |/i| increases. However, the required fine-tuning is still mild e.g. for Mq = 1700 
GeV . 

The coupling between the lightest CP-even Higgs boson and the Z boson is almost 
the same as one in the standard model for small tan 0, however, can decrease to O(10) 
% for moderate tan where the lightest CP-even Higgs is mainly singlet. If the lightest 
CP-even Higgs escapes the LEPII bound, the second lightest CP-even Higgs boson could 
be the boson observed in LHC. The mass of the second lightest CP-even Higgs boson 
depends on k and A, and it can be light in the parameter region favored in our scenario. 
The heaviest CP-even Higgs mass as well as the heaviest CP-odd and charged ones is of 
O(M ). The lightest CP-odd Higgs can also be light. Thus, the Higgs sector has a rich 
structure. 

In our scenario, the higgsino is light compared with three gauginos. In addition, the 
singlino is also light. Both masses of the higgsino and the singlino are of the same order. 
Then, the lightest superparticle is a linear combination between the higgsino and singlino. 
Such a neutralino sector and Higgs sector would lead to several phenomenologically in- 
teresting aspects. 

Note added 

While this work was being completed, we received Ref. [26J, which also considered 
relevant aspects. 
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A Soft SUSY breaking terms 

Here we give explicitly soft SUSY breaking terms induced by the mirage mediation mech- 
anism in the NMSSM. 

In the mirage mediation, the soft parameters at the scale just below Mqvt are given 
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by 



M a (M GUT ) = Mo + ^bagl, 
A ijk (M GUT ) = (a + cj + c fc )M - (7i + 7j + 7 fc ) 
m 2 (M GC/T ) = Cl M 2 - 7j (^) 2 -^M ^, (A.l) 



where 



6 a = -3tr(T a 2 (Adj)) + ^tr(T a 2 (^)), 

i 

1% = 2j29 2 a c^)-lJ2\y^\^ 



2 

jk 



9i = 4j]^C 2 a (^) -^a ijk \y ijk \ 

a jk 



7, = 8n*-p^. (A.2) 
din /i R 

Here, T 2 (Adj) and T 2 (0*) denote Dynkin indices of the adjoint representation and the 
representation of matter fields 0\ We have assumed w„ = T lk iy i j k y* kl to be diagonal. 

Within the framework of the NMSSM, the /3-function coefficients, anomalous dimen- 
sions and other coefficients in the above equations are obtained as 



h = 


-3, b 2 = 1, h = 11, 


1H U = 


\g\ + \gl - 3yt - a 2 , 


lH d = 


2^2 + \al ~ 


Is = 


-2k 2 - 2A 2 , 


IQa = 


+ \al + ^ 2 - (y 2 + yl)&: 


lU a = 


8 n 8 2 or 
3 #3 + gfl - 2j/ t d3a, 


lD a = 


8 o ^9 r\ 2 c 
3^3 + gfl - 2 %^3a, 


lL a = 


3 o -^2 2 r 
2^2 + 2^1 - ^3a, 


lE a = 


2^ 2 - 2 2/ 2 <5 3 a, 



3a • 



(A.3) 
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where 



6h u - 3gj + g\ - 6y 2 aH u Q 3 u§ ~ 2\ 2 a HuHdS , 

2 2 2 2 \ 2 

VH d = 3g 2 + g 1 - Qy b a Hd Q 3 D? i - 2y T a Hd L 3 Ei - 2 A a HuHdS) 
Os = — 2\ 2 aH u H d s — K 2 asss, 

e Qa = ^gl + 3gl + \gl-z(y 2 a HuQzU c + yla HdQ3D c)5z a , 

16 2 16 2 , 

Vu a = y 93 + ~q 9i ~ tClHuQ:iUI 

16 2 4 2 „ 2 
0D O = y#3 + gfl - ^y b a H d Q 3 D-d3a, 

°L a = 3g 2 + g\ - 2y 2 T a HdLsE c5 3a , 

@E a = Ag\ - Ay 2 T a HdL3E cha, (A.4) 



3 11 

in u = 2^2 + y#f - 3y t 2 6 W - A 2 6 A , 

3 4 11 4 o 2t 2t \2t 

= + y#i ~ 3ybbyb ~ Vt Vt ~ X 6a ' 
7s = -2K 2 b K - 2X 2 b x , 
3 11 

iQa = -%3 + 2^2 + ^3^1 ~ ( A* + yb b y b ) 5 3a, 

88 

7c/ = -%s + y#i - tytbytha, 
22 

lD a = ~8#3 + —g\ - 2y 2 b yb 6 3a , 
3 4 11 

7e b = 22gt-2y 2 T b yT 5 3a , (A.5) 



7i a = o^2 + TT^l ~ vlK^a, 



16 q 9 13 o o 9 \ 9 

o w = -y^s -3^2 - y^i +%< +1/6 + A , 

V = -y^-s^-^ + y' + e^ + ^ + A 2 , 

K = -3^ 2 -3^ + 3y 6 2 + 4y 2 + A 2 , 
6 K = 6A 2 + 6/t 2 , 

&a = -3^-(7 2 + 3t/ 2 + 3y 2 + y 2 + 4A 2 + 2 K 2 . (A.6) 

Here, Q a ,U a , D a , L a , and E a denote left-handed quark, right-handed up-sector quark, 
right-handed down-sector quark, left-handed lepton, and right-handed lepton fields, re- 
spectively, and the index a denotes the generation index. We have included effects due to 
Yukawa couplings, yt,yb, and y T , only for the third generations. 
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Figure 1: Constraints and SUSY breaking parameters. We use M = 1200GeV,a = 
2,ch<i = 1,chu = 0, cs = 0, Cq = 0.5, q = 0.5 and tan/3 = 3, A K = — lOOGeV. For detail 
of the shaded regions, see the text. The TeV scale mirage mediation disfavors the region 
\™s I > Mq/V^ ~ 100 GeV (pink shaded). 
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